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ABSTRACT 

We report the detection of a statistically significant flare-like event in the Mg IIA2800A emission line 
of 3C 454.3 during the outburst of autumn 2010. The highest levels of emission line flux recorded over 
the monitoring period (2008 - 2011) coincide with a superluminal jet component traversing through 
the radio core. This finding crucially links the broad-emission line fluctuations to the non-thermal 
continuum emission produced by relativistically moving material in the jet and hence to the presence 
of broad-line region clouds surrounding the radio core. If the radio core were located at several parsecs 
from the central black hole then our results would suggest the presence of broad-line region material 
outside the inner parsec where the canonical broad-line region is envisaged to be located. We briefly 
discuss the implications of broad-emission line material ionized by non-thermal continuum on the 
context of virial black hole mass estimates and gamma-ray production mechanisms. 



1. INTRODUCTION 



It has been long known that the optical emission of 
the flat spectrum r adio quasar 3C 454.3 at z = 0.859 
is highly variable ( |Angione 1968), showi ng structural 



and flux var iability on its parsec scale jet ( Pauliny-Toth 



et al. 1987) and al so correlated variabilit y among dif- 



ferent wave bands (Tornikoski et al. 1994). The recur 



rent flaring behaviour of 3(J 454.3 over the last six years 
(2005-2011, iFuhrmann et al ||2006| iGiommi et aT||2006l 
Raite rTet al.|2008b|a|| Vercellone et al.|2008||Strianl et al. l 
20101 lAckcrmann et al 1120101 IVerccllonc et al.||2011D has 



permitted the assembly of an exquisite multiwavclcngth 



time resolved database (e.g. |Bonning et al.|[2009 Raiteri 



et al. 



varia 



2011 Sasada et al. 12012) allowing to pro be mtraday 



b ility i n the source (e.g [Foschini et al. 2010 Gaur 
et al.]|2012|). Despite the precious database compiled, 



no general consensus about the location of the gamma- 
ray production zone i n 3C 454.3 has been reached so 
far. Some stu dies (e.g. Tavecchio et al. 2010 Poutanen 
& Stern 2010) favor the scenario where the gamma-ray 



emission in 3C 454.3 is generated close to the central 
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black hole (BH) within the canonical broad-line region 
(BLR) , which is locate d within the inner parsec . How- 



ever, other works (e.g. Sikora et al. 2008 Jorstad et al. 



2010[ ) find that the scenario where gamma-rays are pro 



duced far from the central BH, within or downstream of 
the radio core at distances much larger than 1 pc, best 
reproduces the observations of 3C 454.3. 

In this work we explore the variability of the broad- 
emission lines in 3C 454.3 in order to use it as an auxiliary 
piece of information to probe the geometry and physics of 
the innermost regions of 3C 454.3 and to provide evidence 
for the above scenarios of the gamma-ray production. 
Although 3C 454.3 has been monitored extensively at all 



systematically studied dRaiteri et al.||2008a 


Smith et al. 


2009 


Benrtez et al.||2010p. The present wor 


k is the first 



during the Fermi/ L AT era with the largest sample of 
its optical spectra ever compiled. 

2. OBSERVATIONS 

The optical spectra used in this work are taken from 
the Ground-based Observational Support of the Fermi 
Gamma-ray Space Telescope at the University of Ari- 
zona monitoring prograirJ 12 [ details on the o bservational 
setup and reduction process are presented in |Smith et al.| 
(2009). In this work we only consider spectra that have 
been calibrated against the V-band magnitude. The 
spectra have been brought to the rest frame of the source, 
a cosmological correction of the form (1 + z) 3 has been 
performed and no correction for galactic reddening was 
applied. The left panel of Figure [T] displays three optical 
spectra of 3C 454.3 taken at three different periods of 
activity. Here, we present a brief description of our spec- 
tral fitting, which is based on least-squares m inimization 
using the MPFIT package (Markwardt 2009) and closely 
follow s the methodology described in Torrealba et al . 
(2012). Firstly, the featureless continuum is approxi 
mated by a power law function and subtracted from the 
spectrum. Then, the optical Fe II emission is fitted using 



http: / /james.as.arizona.edu/?psmith/Fermi/Website 
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Fig. 1. — Left: Rest frame optical spectra of 3C 454.3 in three intensity states. As can be seen from the top spectrum, the Mg IIA2800A 
emission line is detectable, despite the high levels of optical continuum emission observed. Right: Comparison of the variations seen in the 
Mg IIA2800A profiles (after continuum and Fe II subtraction) from April to November 2010. The observed spectra are shown in solid line 
and a fitted gaussian to the profile are shown in dashed lines. Observing times are color coded as shown in the legend. 

The ratio a/max s is basically the inverse of the S/N ra- 
tio of the line-peak, and the second term Fu ne — F iron i ess 
is introduced to scale the S/N to the fraction of the line 
flux that could be affected by the iron subtraction. The 
third error factor is a result of flux calibration, which 
is approximated by the RMS of flux measurements from 
contiguous observing campaigns (P. Smith, private com- 
munication). Here, we take a very conservative approach 
and approximate the flux calibration error by the stan- 
dard deviation of the flux over the three-year monitoring 
period excluding the observing season showing the high- 
est levels of Mg IIA2800A line flux (> 2a). The sum 
of the first and second error factors is considered here as 
the error of measurements (uF Mg ) and is presented in the 
data points shown in Figures 1 to 4 and listed in Table [lj 
However, the sum of all three sources of error, hereafter 
called total error (a i £ tal ), is used in our statistical anal- 

ysis of variability. The UV-continuum A3000A and Fe II 
fluxes are also listed in Table [TJ 

Figure [2] shows the multiwavelength behaviour of 
3C 454.3 during the Mg IIA2800A monitoring period. 



the template of Vestergaard fc Wilkes| ([2001]) . After the 
local continuum and Fe 11 emission subtraction (see right 
panel of Figure[l]), for all of the spectra, the Mg IIA2800A 
emission line flux was measured by integrating the line 
profile in the range AA2725-2875 A . The spectra m our 
database were taken with different slit sizes (different res- 
olutions) which prevents us to do an immediate analysis 
of the Mg IIA2800A FWHM evolution. Therefore, we 
only measure the flux of Mg II line which is not affected 
by spectral resolution. 

The errors associated with the Mg IIA2800A flux mea- 
surements comprise three different error factors. The 
first error factor is the random error due to disper- 
sion of the spectra and the sign al to noise ratio, esti- 
mated as in |Tresse et al.| ( 1999 ) . The second factor is 
introduced by the subtraction of Fe emission, 
-(flme - F ironless ), where a is the rms of the 
spectrum after iron subtraction, max s is the line-peak in 
the continuum-less spectrum, Fu ne is the line flux, and 
is the fraction of the line flux within the range 



the error 

<JFe = 



F, 



ironless 



of integration where no iron subtraction was performed. 
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taken from the Steward monitoring program. The R- 
band light curve shown in panel (/) is taken from the 
Yale Fermi/SMARTS projeciF] The single dish moni- 
toring data at 37 GHz (A8mmj presented in panel (g) is 
adopted from the Metsah ovi quasar monitoring program 



Fig. 2. — Multiwavelength evolution of 3C 454.3. The vertical stripes show the time when new blobs were ejected from the radio core 
and their widths represent the associated uncertainties. It should be noticed that the highest levels of Mg IIA2800A line flux occurred after 
A8mm flare onset, during an increase in the optical continuum and polarization percentage and before the emergence of a new superluminal 
component from the radio core. 

The gamma-ray light curve from 0.1 to 200 GeV shown 
in panel (a) was built by using data from the Fermi 
Large Area Telescope (LAT) reduced and analyzed with 
the Fermi Science Tools v9r23pl along with the lat- 
est diffuse model files. We have used the user-contributed 
software to generate the spectral models and produce 
the gamma-ray light curve. The Mg IIA2800A flux 
evolution is shown in panel (&). The light curve of 
the Fe II emission is displayed in panel (c). The UV- 
continuum A3000A flux shown in panel (d) has been 
measured from the featureless continuum after removal 
of Fe II template, and the error reported is the RMS 
of the spectrum within 2900-2950 A. The optical lin- 
ear polarization data displayed in panel (e) was also 



now 



( Teraesranta et _&}.._ 
panel of Figure [2] s 
from the radio core 
refer to radio core as 



1998|). The vertical stripes in each 
he ejection tim es of new blobs 



(Jorstad et al 
the unresolved 



2012]. Hereafter, we 
feature in A 7mm 



(43 GHz) VLBI maps associated with the emergence of 
new superluminal components, unless the frequency of 
observation is explicitly stated. 



13 http://www.astro.yale.edu/smarts/glast/tables/3C454.tab 
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Fig. 3.— Flux evolution of Mg IIA2800A emission line (top 
panel) and UV-continuum A3000A emission (bottom panel). For 
each panel, the solid (blue) horizontal line denotes the mean flux 
(fi) observed during the monitoring period, whereas dashed (gray) 
horizontal lines show multiples of cr, where a is the standard de- 
viation of the flux. In this work, we consider a significant flare if 
the levels of emission exceed 2a. Symbol size and color are coded 
according to the color bar displayed, where the larger and lighter 
the symbols, the higher the level of A8mm emission observed. The 
vertical stripes are as in Figure [2] 

3. VARIABILITY OF THE Mg IIA2800A EMISSION LINE. 

The flux evolution of the integrated Mg IIA2800A 
emission line (top panel) and the rest frame UV- 
continuum A3000A (bottom panel) are shown in Figure 
[3] and the data is displayed on Table [T] The significance 
of the fluctuations with respect to the mean flux (//) are 
shown over-plotted as dashed lines. The flux density at 
37 GHz (A8mm) is symbol size and color coded - the 
larger and lighter the symbol, the more intense the mil- 
limeter emission. Variability of the Mg IIA2800A emis- 
sion can be already suggested from visual inspection of 
the light curve and it is confirmed by the large \ 2 values 
obtained with the MPFIT package by fitting a constant 
flux. Under the assumption that the Mg IIA2800A light 
curve can be represented by a constant line, the best fit 
shows a value of x 2 = 291. Hence, taking into account 
the degrees of freedom of the fit (N<j / = 205) , the proba- 
bilities of getting this x 2 value by chance (estimated with 
routine MPCHITEST of MPFIT library) is P < 1 x 1(T 4 . 



Such a small probability allows us to reject the hypothe- 
sis that the emission from the Mg IIA2800A line did not 
vary over the last three years. 

However, no evidence of statistically significant vari- 
ability in the Ha flux of 3C 454.3 was found by |Raiteri 



et al. (2008a I after analyzing 16 NIR spectra taken over a 
six-month period in 2007 (12 June - 03 Nove mber) during 
a faint optical state. Similar to Figure 8 of |Raiteri et al.[ 



(2008a 



we show the response of the Mg IIA2800A line 
flux to the UV-continuum A3000A in Figure [Z] The root 

mean square uncertainty e rms = *J ^ (where are 
the individual errors and N is the number of spectra) 
of the Mg IIA2800A line flux is over-plotted as a dot- 
dashed line. As can be seen, the root mean square of the 
uncertainties in the measurements of the Mg IIA2800A 
line flux is considerably smaller than its standard devi- 
ation (a > e r ms)- F ollowing the same reasoning that 
Raiteri et al. (2008a) used to discard variability in the 
Ha line emission of 3C 454.3, we then find that over 
the period of time considered in this work (2008-2011) 
the Mg IIA2800A line flux variability is well resolved and 
therefore statistically significant. 

We find that during the monitoring period, the spectral 
variation {F max /F min ) of the Mg IIA2800A line is 1.7. 
The peak-to-peak change in the Mg IIA2800A light curve 
is 60% of the mean line flux, while the minimum and 
maximum Mg IIA2800A line flux reached are about 20% 
and 40% of the mean line flux, respectively. Because of 
the sparsely sampled Mg IIA2800A light curve we are 
only able to estimate an upper-limit to its variability 
timescale: r < 1 month. 

Now that we have shown that variability of 
Mg IIA2800A in 3C 454.3 is statistically significant, it 
is instructive to investigate whether there is a causal 
link between major Mg IIA2800A line flux variations and 
variability at other wavelengths, see Figure [2] We have 
also intentionally added a color bar in Fig^f] and [4] to 
show that the highest levels of Mg IIA2800A line flux 
(> 2a), occurred over a short period of time (5500 < 
JD - 2450000 < 5600), coincide with the A8mm fluxes 
above 35 Jy (light red to yellow in color scale). The ver- 
tical stripes in Figure [2] and [3] show the times when new 
blobs were ejected from the radio core, thus suggesting 
that the highest levels of Mg IIA2800A line flux occurred 
after the onset of a mm flare and before a new component 
was seen for the first time leaving the radio core. A mm 
flare presumably start s to rise as the moving dist urbance 
enters the radio core ( Leon-Tavares et al. [2011b ) and by 
definition, the ejection time indicates when such com- 
ponent pa ssed through the middle o f the radio core; see 
Figure 3 in Leon-Tavares et al. (2012). Therefore, the ob- 
served causality of events suggests that the highest levels 
of Mg IIA2800A emission line flux were reached when a 
supcrluminal jet component was traversing through the 
radio core 

How can we be sure that the ionizing continuum was 
produced within the radio core and not by other thermal 
sources like the accretion disk? The high degree of lin- 
ear polarization (see panel e in Figure ul) together with 
the short time scale of variability and the presence of 
a gamma-ray counterpart, suggest that the major UV- 
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Fig. 4.— The response of the Mg IIA2800A emission line (right) and 
emission, where the symbol color-size coding and nomenclature of horiz 
the root mean square uncertainty (e rms ). 

continuum A3000A flare had a non-thermal origin. This 
can be taken as an evidence for a non-thermal source be- 
ing responsible for the ionization of the BLR clouds. In 
further support of this suggestion, in the left panel of Fig- 
ure[4]we can identify a correlation between Mg IIA2800A 
and UV-continuum A3000A (p ~ 0.4) which becomes sig- 
nificant (P < 0.05) when the A8mm emission rises above 
30 Jy. In addition, we also find a significant correla- 
tion between Fe II and UV-continuum A3000A emission 
see right pa nel in Figure HI) in consistency with |Bemtez| 
et al. (2010). The latter Decomes tighter with mcreas- 
mg A8mm emission, thus also supporting the scenario 
in which a complex structure of the BLR is ionized by 
the jet. In order to determine the kinematical proper- 
ties of the BLR clouds illuminated by the inner jet, we 
will perform a detailed analysis of the emission line pro- 
files (studying the FHWM Mg IIA2800A evolution) in a 
forthcoming paper. 

At this point, a discerning reader might ask, why was 
not a Mg IIA2800A flare associated to the component 
ejected in late 2009? Although explaining the physics of 
correlated multiwavelength behaviour in 3C 454.3 is out 
of the scope of this w ork (and we refer the reader to other 
recent studies, e.g. |Bonnoli et al.||2011| |Raiteri et al.| 



2011 Wehrle et al.|2012[ ), we speculate that the detection 
of high levels of Mg IIA2800A emission might be related 
to the amount of energy dissipated by the jet component 
when passing through the radio core. From Figure [2j it 



Fe II emission (right) to levels of UV-continuum A3000A continuum 
ontal lines is the same as in Figure [3] The dot-dashed lines represent 

can be gleaned that the energy dissipation in the events 
of 2009 and 2010 took place downstream and upstream 
of the radio core, respectively. By using the maxima 
of UV-continuum A3000A emission as a proxy for the 
energy dissipated, we see that the event of 2010 released 
twice as much energy than the event from 2009. If the 
amount of energy released at the premises of the radio 
core depends on whether the energy dissipation occurs 
upstream or downstream of the radio core is something 
that deserves further investigation. 

4. SUMMARY AND DISCUSSION 

The optical spectra of 3C 454.3 have been acquired, 
as part of the Ground-based Observational Support of 
the Fermi Gamma-ray Space Telescope at the University 
of Arizona monitoring program, over a period of three 
years (2008-2011) and due to its redshift (z = 0.859) we 
have had access to the middle-UV region of the spectrum, 
allowing us to monitor its Mg IIA2800A broad-emission 
line and adjacent UV-continuum A3000A. We summarize 
our results as follows: 

• We find a statistically significant flare-like event 
in the Mg IIA2800A light curve of 3C 454.3. The 
maximum Mg IIA2800A line flux recorded is above 
40% of the mean flux and the spectral variation 
(F max /F min ) of the Mg IIA2800A line is 1.7. How- 
ever, due to the sparsely sampled Mg IIA2800A 



G 



light curve we are only able to estimate an upper- 
limit to its variability timescale r < 1 month. 

• The highest levels of Mg IIA2800A line flux (> 2er) 
occurred after A8mm flare onset, during an increase 
in the optical polarization percentage, before the 
emergence of a new superluminal component from 
the radio core and within the largest 7-ray flare 
ever seen. This finding crucially links the broad- 
emission line fluctuations to the non-thermal con- 
tinuum emission produced by relativistically mov- 
ing material in the jet and hence to the presence 
of broad-line region clouds surrounding the radio 
core. 

The results presented above indicate that the radio 
core plays a pivotal role in the energy release of 3C 454.3. 
However, there is no general consensus about the true na- 
ture - and location - of the radio core in radio-loud AGN. 
Some studies propose that the radio core is a recollima- 
tion shock and a ge nuine stationary fea ture along the 



parsec scale jet (e.g. Gomez et al. 1995) thought to be 
located at a considerable distance from the black hole. 
On the other hand, recent observational evidence sug- 
gests that for sources with a misaligned jet (i.e. M 87) 



the r adio core mi ght coincide with the BH (Hada et al 
2011| . Recenlty, |Pushkarev et~aL] ( |2012[ ) studied the 
nuclear opacity in the parsec scale jet of a large sam- 
ple of AGN. The authors report that, for 3C 454.3, the 
distance from the radio core at 15 GHz to the apex of 
the jet is about 20 parsecs. Although the radio core at 
43 GHz should be located closer to the BH, the findings 



of Pushkarev et al. (2012) encourage the view that the 
radio core in 3C 454.3 is located well downstream (out- 
side) the canonical BLR. 

Despite the true nature and location of the radio core 
in 3C 454.3 the most important implication of our main 
finding - jet powering broad-emission lines during intense 
outbursts - is the presence of broad-line region material 
surrounding the radio core. If the radio core is embed- 
ded within the canonical BLR, then the ionization of the 
BLR by the jet would make the mirror model to the 
production of gamma-rays as proposed by Ghiscllini & 



Madau ( 1996 ) feasible. Moreover, a likely interaction of 
BLR clouds with the base of the jet might co ntribute 
to the observed levels of high-energy emission (Araudo 
et al.|2010} |Bosch-Ramon et al.|2012[) . Conversely, it the 
radio core is located far from the BH, then broad-line 
region material should be present at distances of parsecs 
from the BH. 

Observational evidence for the presence of BLR mate- 
rial located at parsec scales down the radio core has been 
found by coordinated spectrosc opic and VLBI monitor- 
ing studies. More specifically , Arshakian et al.| (2010) 
and Leon-Tavares e t al. ( |2010 1 found that for the radio 
galaxies 3C 390.3 arid 3C 120 [£] the variable opti- 
cal continuum starts to rise when a new superluminal 
component leaves the radio core seen at 15 GHz and its 
maximum occurs when the component passes through a 
stationary feature located downstream of the radio core. 
Since these two radio galaxies are known to reverberate 

14 http : / /www.metsahovi . f i/~leon/movies/3c3903 . gif 

15 http : //www.metsahovi . f i/-leon/movies/3cl20 .gif 



( Shapovalova et aL]|2010 Grier et aT||2012 ), in the sense 
that the tip broad-emission line responds to changes in 
the optical continuum, the authors conclude that the jet 
can power a significant amount of broad-line emission 
particularly during strong continuum flares in these ob- 
jects. 

The above results suggest the presence of an addi- 
tional component of the BLR, dubbed as outflowing 
BLR, which in effect might be filled with BLR mate- 
rial dragged by the relativistic jet as it propagates down- 
stream of the BH or perhaps could be a sub-relativistic 
outflow arising from an accretion-disk wind. The notion 
of a dynamic and extended BLR has been previously pro 



posed in several works (e.g. Popovic et aLpOOl 
|Shlosma n|2006 ) . Although these studies nave 



Elitzur fc 
ocused on 

non-blazar sources, the presence of an outflowing BLR 
has already been su ggested for so urces like 3C 273 flPal 



tani fc Turler]2003| ) and 3C 454.3 ( |Finke k. Dermer|2i 
where it could serve as an alternative source of seed pho 
tons for the Inverse Compt on scattering as proposed by 



Leon-Tavares et al. (2011b). 

Regardless of whether the radio core is close (1 pc < ) 
or far (^1 pc) from the BH, the fact that broad-emission 
lines respond to changes of the non-thermal continuum 
prevents us to use the single epoch virial BH mass es- 
timates in this source because the latter assumes: (i) a 
single localized ionization source (i.e. accretion disk) and 
(ii) virial equilibrium of the BLR clouds. The latter as- 
sumptions cannot be fulfilled during episodes of strong 
flaring activity in 3C 454.3, hence the ionization of BLR 
clouds by non-thermal emission might introduce uncer- 
tainties to the BH mass estimates derived by assuming 
virial equilibrium of the BLR. Alternative scaling rela- 
tions to weigh the BH in str ongly beamed sources are dis - 



cussed and implemented in |Le6n-Tavares et al. (2011a 
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TABLE 1 

Measurements 



JD-2450000 


Mg IIA2800A 




Fe II 




UV-continuum A3000A 






flux (10 -13 erg cm -2 s _1 ) 


error flux (10 


-13 er g cm -2 g — 1 ) 


error 


flux(10" 14 erg cm- 2 s" 1 A" 1 ) 


error 


4743.8350 


1.412 


0.019 


2.418 


0.478 


1.759 


0.005 


4745.7251 


1.261 


0.016 


1.962 


0.476 


1.698 


0.005 


4746.7104 


1.248 


0.018 


2.094 


0.477 


1.648 


0.005 


4747.8491 


1.199 


0.016 


2.333 


0.474 


1.850 


0.005 


4748.7373 


1.336 


0.028 


2.164 


0.479 


2.159 


0.009 


4767.6724 


1.322 


0.016 


2.389 


0.473 


0.864 


0.005 


4768.6597 


1.357 


0.018 


2.416 


0.473 


0.971 


0.006 


4769.8086 


1.321 


0.014 


2.308 


0.474 


0.923 


0.004 


4770.7051 


1.320 


0.016 


2.101 


0.475 


0.821 


0.004 


4772.6831 


1.327 


0.014 


2.020 


0.474 


0.788 


0.004 



Note. — Tablc^is published in its entirely in the electronic edition 
of the Astrophysical Journal Letters. A portion is shown here for 
guidance regarding its form and content. 
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